Hamstrings are stretched more and faster during accelerative running compared to speed-matched constant speed running.

Background: Hamstring strain injuries are associated with significant time away from sport and high reinjury rates. Recent evidence suggests that hamstring injuries often occur during accelerative running, but investigations of hamstring mechanics have primarily examined constant speed running on a treadmill. To help fill this gap in knowledge, this study compares hamstring lengths and lengthening velocities between accelerative running and constant speed overground running. Methods: We recorded 2 synchronized videos of 10 participants (5 female, 5 male) during 6 accelerative running trials and 6 constant speed running trials. We used OpenCap (a markerless motion capture system) to estimate body segment kinematics for each trial and a 3-dimensional musculoskeletal model to compute peak length and step-average lengthening velocity of the biceps femoris (long head) muscle-tendon unit. To compare running conditions, we used linear mixed regression models with running speed (normalized by the subject-specific maximum) as the independent variable. Results: At running speeds below 75% of top speed accelerative running resulted in greater peak lengths than constant speed running. For example, the peak hamstring muscle-tendon length when a person accelerated from running at only 50% of top speed was equivalent to running at a constant 88% of top speed. Lengthening velocities were greater during accelerative running at all running speeds. Differences in hip flexion kinematics primarily drove the greater peak muscle-tendon lengths and lengthening velocities observed in accelerative running. Conclusion: Hamstrings are subjected to longer muscle-tendon lengths and faster lengthening velocities in accelerative running compared to constant speed running. This provides a biomechanical explanation for the observation that hamstring strain injuries often occur during acceleration. Our results suggest coaches who monitor exposure to high-risk circumstances (long lengths, fast lengthening velocities) should consider the accelerative nature of running in addition to running speed.

OpenCap: Human movement dynamics from smartphone videos.

Measures of human movement dynamics can predict outcomes like injury risk or musculoskeletal disease progression. However, these measures are rarely quantified in large-scale research studies or clinical practice due to the prohibitive cost, time, and expertise required. Here we present and validate OpenCap, an open-source platform for computing both the kinematics (i.e., motion) and dynamics (i.e., forces) of human movement using videos captured from two or more smartphones. OpenCap leverages pose estimation algorithms to identify body landmarks from videos; deep learning and biomechanical models to estimate three-dimensional kinematics; and physics-based simulations to estimate muscle activations and musculoskeletal dynamics. OpenCap's web application enables users to collect synchronous videos and visualize movement data that is automatically processed in the cloud, thereby eliminating the need for specialized hardware, software, and expertise. We show that OpenCap accurately predicts dynamic measures, like muscle activations, joint loads, and joint moments, which can be used to screen for disease risk, evaluate intervention efficacy, assess between-group movement differences, and inform rehabilitation decisions. Additionally, we demonstrate OpenCap's practical utility through a 100-subject field study, where a clinician using OpenCap estimated musculoskeletal dynamics 25 times faster than a laboratory-based approach at less than 1% of the cost. By democratizing access to human movement analysis, OpenCap can accelerate the incorporation of biomechanical metrics into large-scale research studies, clinical trials, and clinical practice.

Smartphone videos of the sit-to-stand test predict osteoarthritis and health outcomes in a nationwide study.

Physical function decline due to aging or disease can be assessed with quantitative motion analysis, but this currently requires expensive laboratory equipment. We introduce a self-guided quantitative motion analysis of the widely used five-repetition sit-to-stand test using a smartphone. Across 35 US states, 405 participants recorded a video performing the test in their homes. We found that the quantitative movement parameters extracted from the smartphone videos were related to a diagnosis of osteoarthritis, physical and mental health, body mass index, age, and ethnicity and race. Our findings demonstrate that at-home movement analysis goes beyond established clinical metrics to provide objective and inexpensive digital outcome metrics for nationwide studies.

Modeling toes contributes to realistic stance knee mechanics in three-dimensional predictive simulations of walking.

Physics-based predictive simulations have been shown to capture many salient features of human walking. Yet they often fail to produce realistic stance knee and ankle mechanics. While the influence of the performance criterion on the predicted walking pattern has been previously studied, the influence of musculoskeletal mechanics has been less explored. Here, we investigated the influence of two mechanical assumptions on the predicted walking pattern: the complexity of the foot model and the stiffness of the Achilles tendon. We found, through three-dimensional muscle-driven predictive simulations of walking, that modeling the toes, and thus using two-segment instead of single-segment foot models, contributed to robustly eliciting physiological stance knee flexion angles, knee extension torques, and knee extensor activity. Modeling toes also slightly decreased the first vertical ground reaction force peak, increasing its agreement with experimental data, and improved stance ankle kinetics. It nevertheless slightly worsened predictions of ankle kinematics. Decreasing Achilles tendon stiffness improved the realism of ankle kinematics, but there remain large discrepancies with experimental data. Overall, this simulation study shows that not only the performance criterion but also mechanical assumptions affect predictive simulations of walking. Improving the realism of predictive simulations is required for their application in clinical contexts. Here, we suggest that using more complex foot models might contribute to such realism.

Predictive simulations of running gait reveal a critical dynamic role for the tail in bipedal dinosaur locomotion.

Locomotion has influenced the ecology, evolution, and extinction of species throughout history, yet studying locomotion in the fossil record is challenging. Computational biomechanics can provide novel insight by mechanistically relating observed anatomy to whole-animal function and behavior. Here, we leverage optimal control methods to generate the first fully predictive, three-dimensional, muscle-driven simulations of locomotion in an extinct terrestrial vertebrate, the bipedal non-avian theropod dinosaur Coelophysis. Unexpectedly, our simulations involved pronounced lateroflexion movements of the tail. Rather than just being a static counterbalance, simulations indicate that the tail played a crucial dynamic role, with lateroflexion acting as a passive, physics-based mechanism for regulating angular momentum and improving locomotor economy, analogous to the swinging arms of humans. We infer this mechanism to have existed in many other bipedal non-avian dinosaurs as well, and our methodology provides new avenues for exploring the functional diversity of dinosaur tails in the future.

Computational modelling of muscle fibre operating ranges in the hindlimb of a small ground bird (Eudromia elegans), with implications for modelling locomotion in extinct species.

The arrangement and physiology of muscle fibres can strongly influence musculoskeletal function and whole-organismal performance. However, experimental investigation of muscle function during in vivo activity is typically limited to relatively few muscles in a given system. Computational models and simulations of the musculoskeletal system can partly overcome these limitations, by exploring the dynamics of muscles, tendons and other tissues in a robust and quantitative fashion. Here, a high-fidelity, 26-degree-of-freedom musculoskeletal model was developed of the hindlimb of a small ground bird, the elegant-crested tinamou (Eudromia elegans, ~550 g), including all the major muscles of the limb (36 actuators per leg). The model was integrated with biplanar fluoroscopy (XROMM) and forceplate data for walking and running, where dynamic optimization was used to estimate muscle excitations and fibre length changes throughout both gaits. Following this, a series of static simulations over the total range of physiological limb postures were performed, to circumscribe the bounds of possible variation in fibre length. During gait, fibre lengths for all muscles remained between 0.5 to 1.21 times optimal fibre length, but operated mostly on the ascending limb and plateau of the active force-length curve, a result that parallels previous experimental findings for birds, humans and other species. However, the ranges of fibre length varied considerably among individual muscles, especially when considered across the total possible range of joint excursion. Net length change of muscle-tendon units was mostly less than optimal fibre length, sometimes markedly so, suggesting that approaches that use muscle-tendon length change to estimate optimal fibre length in extinct species are likely underestimating this important parameter for many muscles. The results of this study clarify and broaden understanding of muscle function in extant animals, and can help refine approaches used to study extinct species.

A Dynamic Optimization Approach for Solving Spine Kinematics While Calibrating Subject-Specific Mechanical Properties.

This study aims to propose a new optimization framework for solving spine kinematics based on skin-mounted markers and estimate subject-specific mechanical properties of the intervertebral joints. The approach enforces dynamic consistency in the entire skeletal system over the entire time-trajectory while personalizing spinal stiffness. 3D reflective markers mounted on ten vertebrae during spine motions were measured in ten healthy volunteers. Biplanar X-rays were taken during neutral stance of the subjects wearing the markers. Calculated spine kinematics were compared to those calculated using inverse kinematics (IK) and IK with imposed generic kinematic constraints. Calculated spine kinematics compared well with standing X-rays, with average root mean square differences of the vertebral body center positions below 10.1 mm and below [Formula: see text] for joint orientation angles. For flexion/extension and lateral bending, the lumbar rotation distribution patterns, as well as the ranges of rotations matched in vivo literature data. The approach outperforms state-of-art IK and IK with constraints methods. Calculated ratios reflect reduced spinal stiffness in low-resistance zone and increased stiffness in high-resistance zone. The patterns of calibrated stiffness were consistent with previously reported experimentally determined patterns. This approach will further our insight into spinal mechanics by increasing the physiological representativeness of spinal motion simulations.

Perspective on musculoskeletal modelling and predictive simulations of human movement to assess the neuromechanics of gait.

Locomotion results from complex interactions between the central nervous system and the musculoskeletal system with its many degrees of freedom and muscles. Gaining insight into how the properties of each subsystem shape human gait is challenging as experimental methods to manipulate and assess isolated subsystems are limited. Simulations that predict movement patterns based on a mathematical model of the neuro-musculoskeletal system without relying on experimental data can reveal principles of locomotion by elucidating cause-effect relationships. New computational approaches have enabled the use of such predictive simulations with complex neuro-musculoskeletal models. Here, we review recent advances in predictive simulations of human movement and how those simulations have been used to deepen our knowledge about the neuromechanics of gait. In addition, we give a perspective on challenges towards using predictive simulations to gain new fundamental insight into motor control of gait, and to help design personalized treatments in patients with neurological disorders and assistive devices that improve gait performance. Such applications will require more detailed neuro-musculoskeletal models and simulation approaches that take uncertainty into account, tools to efficiently personalize those models, and validation studies to demonstrate the ability of simulations to predict gait in novel circumstances.

OpenSim Moco: Musculoskeletal optimal control.

Musculoskeletal simulations are used in many different applications, ranging from the design of wearable robots that interact with humans to the analysis of patients with impaired movement. Here, we introduce OpenSim Moco, a software toolkit for optimizing the motion and control of musculoskeletal models built in the OpenSim modeling and simulation package. OpenSim Moco uses the direct collocation method, which is often faster and can handle more diverse problems than other methods for musculoskeletal simulation. Moco frees researchers from implementing direct collocation themselves-which typically requires extensive technical expertise-and allows them to focus on their scientific questions. The software can handle a wide range of problems that interest biomechanists, including motion tracking, motion prediction, parameter optimization, model fitting, electromyography-driven simulation, and device design. Moco is the first musculoskeletal direct collocation tool to handle kinematic constraints, which enable modeling of kinematic loops (e.g., cycling models) and complex anatomy (e.g., patellar motion). To show the abilities of Moco, we first solved for muscle activity that produced an observed walking motion while minimizing squared muscle excitations and knee joint loading. Next, we predicted how muscle weakness may cause deviations from a normal walking motion. Lastly, we predicted a squat-to-stand motion and optimized the stiffness of an assistive device placed at the knee. We designed Moco to be easy to use, customizable, and extensible, thereby accelerating the use of simulations to understand the movement of humans and other animals.

Physics-Based Simulations to Predict the Differential Effects of Motor Control and Musculoskeletal Deficits on Gait Dysfunction in Cerebral Palsy: A Retrospective Case Study.

Physics-based simulations of walking have the theoretical potential to support clinical decision-making by predicting the functional outcome of treatments in terms of walking performance. Yet before using such simulations in clinical practice, their ability to identify the main treatment targets in specific patients needs to be demonstrated. In this study, we generated predictive simulations of walking with a medical imaging based neuro-musculoskeletal model of a child with cerebral palsy presenting crouch gait. We explored the influence of altered muscle-tendon properties, reduced neuromuscular control complexity, and spasticity on gait dysfunction in terms of joint kinematics, kinetics, muscle activity, and metabolic cost of transport. We modeled altered muscle-tendon properties by personalizing Hill-type muscle-tendon parameters based on data collected during functional movements, simpler neuromuscular control by reducing the number of independent muscle synergies, and spasticity through delayed muscle activity feedback from muscle force and force rate. Our simulations revealed that, in the presence of aberrant musculoskeletal geometries, altered muscle-tendon properties rather than reduced neuromuscular control complexity and spasticity were the primary cause of the crouch gait pattern observed for this child, which is in agreement with the clinical examination. These results suggest that muscle-tendon properties should be the primary target of interventions aiming to restore an upright gait pattern for this child. This suggestion is in line with the gait analysis following muscle-tendon property and bone deformity corrections. Future work should extend this single case analysis to more patients in order to validate the ability of our physics-based simulations to capture the gait patterns of individual patients pre- and post-treatment. Such validation would open the door for identifying targeted treatment strategies with the aim of designing optimized interventions for neuro-musculoskeletal disorders.

Algorithmic differentiation improves the computational efficiency of OpenSim-based trajectory optimization of human movement.

Algorithmic differentiation (AD) is an alternative to finite differences (FD) for evaluating function derivatives. The primary aim of this study was to demonstrate the computational benefits of using AD instead of FD in OpenSim-based trajectory optimization of human movement. The secondary aim was to evaluate computational choices including different AD tools, different linear solvers, and the use of first- or second-order derivatives. First, we enabled the use of AD in OpenSim through a custom source code transformation tool and through the operator overloading tool ADOL-C. Second, we developed an interface between OpenSim and CasADi to solve trajectory optimization problems. Third, we evaluated computational choices through simulations of perturbed balance, two-dimensional predictive simulations of walking, and three-dimensional tracking simulations of walking. We performed all simulations using direct collocation and implicit differential equations. Using AD through our custom tool was between 1.8 ± 0.1 and 17.8 ± 4.9 times faster than using FD, and between 3.6 ± 0.3 and 12.3 ± 1.3 times faster than using AD through ADOL-C. The linear solver efficiency was problem-dependent and no solver was consistently more efficient. Using second-order derivatives was more efficient for balance simulations but less efficient for walking simulations. The walking simulations were physiologically realistic. These results highlight how the use of AD drastically decreases computational time of trajectory optimization problems as compared to more common FD. Overall, combining AD with direct collocation and implicit differential equations decreases the computational burden of trajectory optimization of human movement, which will facilitate their use for biomechanical applications requiring the use of detailed models of the musculoskeletal system.

SimCP: A Simulation Platform to Predict Gait Performance Following Orthopedic Intervention in Children With Cerebral Palsy.

Gait deficits in cerebral palsy (CP) are often treated with a single-event multi-level surgery (SEMLS). Selecting the treatment options (combination of bony and soft tissue corrections) for a specific patient is a complex endeavor and very often treatment outcome is not satisfying. A deterioration in 22.8% of the parameters describing gait performance has been reported and there is need for additional surgery in 11% of the patients. Computational simulations based on musculoskeletal models that allow clinicians to test the effects of different treatment options before surgery have the potential to drastically improve treatment outcome. However, to date, no such simulation and modeling method is available. Two important challenges are the development of methods to include patient-specific neuromechanical impairments into the models and to simulate the effect of different surgical procedures on post-operative gait performance. Therefore, we developed the SimCP framework that allows the evaluation of the effect of different simulated surgeries on gait performance of a specific patient and includes a graphical user interface (GUI) that enables performing virtual surgery on the models. We demonstrated the potential of our framework for two case studies. Models reflecting the patient-specific musculoskeletal geometry and muscle properties are generated based solely on data collected before the treatment. The patient's motor control is described based on muscle synergies derived from pre-operative EMG. The GUI is then used to modify the musculoskeletal properties according to the surgical plan. Since SEMLS does not affect motor control, the same motor control model is used to define gait performance pre- and post-operative. We use the capability gap (CG), i.e., the difference between the joint moments needed to perform healthy walking and the joint moments the personalized model can generate, to quantify gait performance. In both cases, the CG was smaller post- then pre-operative and this was in accordance with the measured change in gait kinematics after treatment.

Rapid predictive simulations with complex musculoskeletal models suggest that diverse healthy and pathological human gaits can emerge from similar control strategies.

Physics-based predictive simulations of human movement have the potential to support personalized medicine, but large computational costs and difficulties to model control strategies have limited their use. We have developed a computationally efficient optimal control framework to predict human gaits based on optimization of a performance criterion without relying on experimental data. The framework generates three-dimensional muscle-driven simulations in 36 min on average-more than 20 times faster than existing simulations-by using direct collocation, implicit differential equations and algorithmic differentiation. Using this framework, we identified a multi-objective performance criterion combining energy and effort considerations that produces physiologically realistic walking gaits. The same criterion also predicted the walk-to-run transition and clinical gait deficiencies caused by muscle weakness and prosthesis use, suggesting that diverse healthy and pathological gaits can emerge from the same control strategy. The ability to predict the mechanics and energetics of a broad range of gaits with complex three-dimensional musculoskeletal models will allow testing novel hypotheses about gait control and hasten the development of optimal treatments for neuro-musculoskeletal disorders.

Subject-Exoskeleton Contact Model Calibration Leads to Accurate Interaction Force Predictions.

Knowledge of human-exoskeleton interaction forces is crucial to assess user comfort and effectiveness of the interaction. The subject-exoskeleton collaborative movement and its interaction forces can be predicted in silico using computational modeling techniques. We developed an optimal control framework that consisted of three phases. First, the foot-ground (Phase A) and the subject-exoskeleton (Phase B) contact models were calibrated using three experimental sit-to-stand trials. Then, the collaborative movement and the subject-exoskeleton interaction forces, of six different sit-to-stand trials were predicted (Phase C). The results show that the contact models were able to reproduce experimental kinematics of calibration trials (mean root mean square differences - RMSD - coordinates ≤ 1.1° and velocities ≤ 6.8°/s), ground reaction forces (mean RMSD≤ 22.9 N), as well as the interaction forces at the pelvis, thigh, and shank (mean RMSD ≤ 5.4 N). Phase C could predict the collaborative movements of prediction trials (mean RMSD coordinates ≤ 3.5° and velocities ≤ 15.0°/s), and their subject-exoskeleton interaction forces (mean RMSD ≤ 13.1° N). In conclusion, this optimal control framework could be used while designing exoskeletons to have in silico knowledge of new optimal movements and their interaction forces.

A spasticity model based on feedback from muscle force explains muscle activity during passive stretches and gait in children with cerebral palsy.

Muscle spasticity is characterized by exaggerated stretch reflexes and affects about 85% of the children with cerebral palsy. However, the mechanisms underlying spasticity and its influence on gait are not well understood. Here, we first aimed to model the response of spastic hamstrings and gastrocnemii in children with cerebral palsy to fast passive stretches. Then, we evaluated how the model applied to gait. We developed three models based on exaggerated proprioceptive feedback. The first model relied on feedback from muscle fiber length and velocity (velocity-related model), the second model relied on feedback from muscle fiber length, velocity, and acceleration (acceleration-related model), and the third model relied on feedback from muscle force and its first time derivative (force-related model). The force-related model better reproduced measured hamstrings and gastrocnemii activity during fast passive stretches (coefficients of determination (R2): 0.73 ± 0.10 and 0.60 ± 0.13, respectively, and root mean square errors (RMSE): 0.034 ± 0.031 and 0.009 ± 0.007, respectively) than the velocity-related model (R2: 0.46 ± 0.15 and 0.07 ± 0.13, and RMSE: 0.053 ± 0.051 and 0.015 ± 0.009), and the acceleration-related model (R2: 0.47 ± 0.15 and 0.09 ± 0.14, and RMSE: 0.052 ± 0.050 and 0.015 ± 0.008). Additionally, the force-related model predicted hamstrings and gastrocnemii activity that better correlated with measured activity during gait (cross correlations: 0.82 ± 0.09 and 0.85 ± 0.06, respectively) than the activity predicted by the velocity-related model (cross correlations: 0.49 ± 0.17 and 0.71 ± 0.22) and the acceleration-related model (cross correlations: 0.51 ± 0.16 and 0.67 ± 0.20). Our results therefore suggest that force encoding in muscle spindles in combination with altered feedback gains and thresholds underlie activity of spastic muscles during passive stretches and gait. Our model of spasticity opens new perspectives for studying movement impairments due to spasticity through simulation.

The influence of maximum isometric muscle force scaling on estimated muscle forces from musculoskeletal models of children with cerebral palsy.

BACKGROUND: Musculoskeletal models do not include patient-specific muscle forces but rely on a scaled generic model, with muscle forces left unscaled in most cases. However, to use musculoskeletal simulations to inform clinical decision-making in children with cerebral palsy (CP), inclusion of subject-specific muscle forces is of utmost importance in order to represent each child's compensation mechanisms introduced through muscle weakness. RESEARCH AIM: The aims of this study were to (i) evaluate if maximum isometric muscle forces (MIMF) in musculoskeletal models of children with CP can be scaled based on strength measurements obtained with a hand-held-dynamometer (HHD), (ii) evaluate the impact of the HHD based scaling approach and previously published MIMF scaling methods on computed muscle forces during gait, and (iii) compare maximum muscle forces during gait between CP and typically developing (TD) children. METHODS: Strength and motion capture data of six CP and motion capture data of six TD children were collected. The HHD measurements to obtain hip, knee and ankle muscle strength were simulated in OpenSim and used to modify MIMF of the 2392-OpenSim model. These muscle forces were compared to the MIMF scaled on the child's body mass and a scaling approach, which included the body mass and muscle-tendon lengths. OpenSim was used to calculate peak muscle forces during gait. RESULTS: Ankle muscle strength was insufficient to reproduce joint moments during walking when MIMF were scaled based on HHD. During gait, peak hip and knee extensor muscle forces were higher and peak ankle dorsi-flexor forces were lower in CP compared to TD participants. SIGNIFICANCE: HHD measurements can be used to scale MIMF for the hip and knee muscle groups but underestimate the force capacity of the ankle muscle groups during walking. Muscle-tendon-length and mass based scaling methods affected muscle activations but had little influence on peak muscle forces during gait.

OpenSim Versus Human Body Model: A Comparison Study for the Lower Limbs During Gait.

Musculoskeletal modeling and simulations have become popular tools for analyzing human movements. However, end users are often not aware of underlying modeling and computational assumptions. This study investigates how these assumptions affect biomechanical gait analysis outcomes performed with Human Body Model and the OpenSim gait2392 model. The authors compared joint kinematics, kinetics, and muscle forces resulting from processing data from 7 healthy adults with both models. Although outcome variables had similar patterns, there were statistically significant differences in joint kinematics (maximal difference: 9.8° [1.5°] in sagittal plane hip rotation), kinetics (maximal difference: 0.36 [0.10] N·m/kg in sagittal plane hip moment), and muscle forces (maximal difference: 8.51 [1.80] N/kg for psoas). These differences might be explained by differences in hip and knee joint center locations up to 2.4 (0.5) and 1.9 (0.2) cm in the posteroanterior and inferosuperior directions, respectively, and by the offset in pelvic reference frames of about 10° around the mediolateral axis. The choice of model may not influence the conclusions in clinical settings, where the focus is on interpreting deviations from the reference data, but it will affect the conclusions of mechanical analyses in which the goal is to obtain accurate estimates of kinematics and loading.

EMG-Driven Optimal Estimation of Subject-SPECIFIC Hill Model Muscle-Tendon Parameters of the Knee Joint Actuators.

OBJECTIVE: the purpose of this paper is to propose an optimal control problem formulation to estimate subject-specific Hill model muscle-tendon (MT-) parameters of the knee joint actuators by optimizing the fit between experimental and model-based knee moments. Additionally, this paper aims at determining which sets of functional motions contain the necessary information to identify the MT-parameters. METHODS: the optimal control and parameter estimation problem underlying the MT-parameter estimation is solved for subject-specific MT-parameters via direct collocation using an electromyography-driven musculoskeletal model. The sets of motions containing sufficient information to identify the MT-parameters are determined by evaluating knee moments simulated based on subject-specific MT-parameters against experimental moments. RESULTS: the MT-parameter estimation problem was solved in about 30 CPU minutes. MT-parameters could be identified from only seven of the 62 investigated sets of motions, underlining the importance of the experimental protocol. Using subject-specific MT-parameters instead of more common linearly scaled MT-parameters improved the fit between inverse dynamics moments and simulated moments by about 30% in terms of the coefficient of determination (from [Formula: see text] to [Formula: see text]) and by about 26% in terms of the root mean square error (from [Formula: see text] to [Formula: see text] ). In particular, subject-specific MT-parameters of the knee flexors were very different from linearly scaled MT-parameters. CONCLUSION: we introduced a computationally efficient optimal control problem formulation and provided guidelines for designing an experimental protocol to estimate subject-specific MT-parameters improving the accuracy of motion simulations. SIGNIFICANCE: the proposed formulation opens new perspectives for subject-specific musculoskeletal modeling, which might be beneficial for simulating and understanding pathological motions.

Preliminary Testing of a Compact, Bone-Attached Robot for Otologic Surgery.

Otologic surgery often involves a mastoidectomy procedure, in which part of the temporal bone is milled away in order to visualize critical structures embedded in the bone and safely access the middle and inner ear. We propose to automate this portion of the surgery using a compact, bone-attached milling robot. A high level of accuracy is required to avoid damage to vital anatomy along the surgical path, most notably the facial nerve, making this procedure well-suited for robotic intervention. In this study, several of the design considerations are discussed and a robot design and prototype are presented. The prototype is a 4 degrees-of-freedom robot similar to a four-axis milling machine that mounts to the patient's skull. A positioning frame, containing fiducial markers and attachment points for the robot, is rigidly attached to the skull of the patient, and a CT scan is acquired. The target bone volume is manually segmented in the CT by the surgeon and automatically converted to a milling path and robot trajectory. The robot is then attached to the positioning frame and is used to drill the desired volume. The accuracy of the entire system (image processing, planning, robot) was evaluated at several critical locations within or near the target bone volume with a mean free space accuracy result of 0.50 mm or less at all points. A milling test in a phantom material was then performed to evaluate the surgical workflow. The resulting milled volume did not violate any critical structures.